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DISCUSISION | 


VEN Te Cnow. Following the computation procedure p proposed in the 


paper, the variability index can be expressed ~s the formula: * 


_ duration curve; g is the mean of y; y — @ is the deviation from the mean; and 
nis the number of selected discharges, or n° = 10i in the given procedure. * ‘hee a 
Ss, operation, it is seldom convenient ‘to take the deviations from the 

actual mean, since ‘such deviations usually involve decimals which are cumber- 
some to handle when squared. For greatest convenience (especially in machine 
Ee form of Eq. 1 is obtained by ut utilizing the mathematical — 


which y is is the ae the selected at 10% interval of 


principle that the sum of the squares of deviation from the mean n is equal to othe 
sum of the squares of deviation from zero, less the product of of the total and the a 
mean. . M Modified slightly, this this principle i is expressed as 


Eq. 2 is known as she Gaussian feeanin, which i is used to estimate the standard 
_ deviation of in statistics. In Eq. 2 yi is the mean of the ‘squares 


The ‘numerical method of computing the variability index was 
presented for the purpose of eliminating the effects introduced by variations - 
in personal judgment. However, as the method used the values of discharges: 
a _ read from the duration curve as s plotted, it did not eliminate the personal factor | 
introduced in plotting: the curve. © A ‘complete elimination of the influence of 
personal judgment is possible ‘only by replacing the ‘ ‘eye fitting” by a mathe- 
matical process— —that i is, the values of : 9 and y contained i in Eq. 2 should be 
computed directly | from the Taw data. Such data can be grouped by the 
method of statistical classification to reduce the number of data handled in 
computation. — if necessary, the c computed standard deviation can be corrected . 
by applying the Sheppard correction. 


_ Since the mathematical part « of the paper | has been omitted, ‘the method 4 ay 


- _Notge.—This paper by E. W. Lane and Kai Lei was published i in September, 1949, Proceedings. ” The 4 
numbering of footnotes, tables, equations, and illustrations in in this Separate is a continuation of the: con- ve 
secutive numbering used in the original paper. 
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CHOW ON: STREAM FLOW 


the mean annual rt runo cannot be checked readily. However, the st stipulation — 


(under the heading, “Use of the Variability Index”) that— 


the ratio of the discharge exceeded 15. 87% of the time to 


discharge exceeded 50% of the time was equal to the variability index 


selected” 


_—would make t the variability index; greater than unity: for most streams, because 
the discharge that is exceeded 15. 877% of the time is usually greater than ‘the 


with the of the variability index listed it in 1 Table 1. 

‘The knowledge of hydrologic statistics has been developed and extended — _ 
"greatly in in recent years. investigators'”: ‘2 have found that the distribu- 
of most hydrologic phenomena is skew. Probability plottings of the 
Hazen type were also found unsatisfactory for such frequencies. The computa- 


tion of asymmetric probabilities i is now possible in the field of hydrology; and eae 


development of a special probability paper!6.23 ‘opened the way fora straight- 
tine plotting of the duration curve. "Therefore the study of stream flow 
would be advanced if a suitable asymmetric probability plotting could be . 


| adapted to it. The values of discharge should be plotted on a linear scale, so. 
_ that the variability index would be equal to the standard deviation, in’ in a 


hydrologists like to use the coefficient of variation, C., instead of “4 
standard deviation, ¢, for the study of variability. The coefficient of 
variation is equal to the standard deviation divided by the mean, or 


4 


"effects of the difference from the mean. This elin elimination is desirable, 
2 stream discharge i in more arid 1 regions where the annual rainfall tends to be ; 


smaller as well as more variable. 
- ie For any st statistical study, the quantity y of the available data i is as as ‘important ; 7 
the quality. As long ai as the q quality i is kept homogeneous, the longer record, 


or larger number of data, would make possible more accurate estimates of the : sé 
; frequencies of stream flow. W. D. Potter,* of the Soil Conservation Service, _ — 


United States Department of Agriculture, has made a probability study 
rainfall and runoff. The study WAS confined to annual | maximum \ values; but 
‘it is worthy to note that the frequency of distribution is a function of years of 
record, Ww hich is a measure of the number of data used . Since the variability 
ae index i is a measure of the slope of the duration curve (a frequency curve), it 


should be some function of the length of record. Therefore, the variability 
_ index would be of more statistical value if _ two were compared on the basis 


2 “Floods Estimated by Probability Method,” "by E. J. News-Record, Sune 14, 

«16 “Straight Line Plotting. of Ghow Data,” by. R. D. Transactions, ASCE, Vol. 91, 

‘‘A Simple Method of Estimating Flood Frequency,” by Ralph]W. Powell, Civil cee, Febru- 


“Simplification of the Gumbel Method for | Computing | Probability Curves,” by w. D. 
No. SCS-TP-78, Soil Conservation Service, U.S. D. A., Washington, D. C., May, 1949. 
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WOOD ON STREAM FLOW 
Dana M. Woon,* M. ASCE.—The selection of an analytical 
; _ _ depends upon the use to which the results are to be put and the time available a 
the study. = Cost i is a also a factor i in any engineering office. Two ‘sugges- 


tions con concerning the duration curve ma may y be of value to those desiring to extend ; 


(1) — of 1 using : a scale of ‘ ‘discharge i in cubic feet 1 per second per ‘square vy 
_ nile,” substituting “ratio to the average flow” for the period of record under | 
eo will result i in relative variability curves of flow. The values 0 of . 7 


departures : at 10% intervals can then be obtained ‘similarly from this ¢ curve. 


runoff i in any comparisons. stream having ¢ a flow of 0. 5 ct cu per sec 
7 _ ~per sq ‘mile can be « compared with one having a flow 2.5 cu ft per sec per sq mile 
a with respect to the slopes of the duration curves—that i is, the magnitude of ‘the | 


departures from the average, regardless of the size and other characteristics of 
7 the drainage area. The writer would like to know if this approach might in 

any way simplify the logarithmic calculations, 

; _ (2) The second point of importance is that the same — of years should ; 


be used when ‘comparing different records in the same general locality. Itis 
; well known that the ‘magnitude of the runoff varies in cycles of wet and dry a 
; years, although little enough is known about the amplitude and magnitude of 
7 the curve of variation. — However, in the Tennessee basin and numerous =a 
= the years from 1930 to 1948 have been generally i in the “dry” part of a 
eycle. The previous nineteen years, from 1910 to 1929, were generally i ina 
flat part of the cycle having normal annual variations. It is quite obvious | 
* - that any analysis could lead to erronous conclusions if records in the two sepa- _ 
7 rate: periods were compared. This point is well illustrated by the doubt | 
regarding the analysis of in different. periods of years when attempting 4 
_ determine the frequency of occurrence of floods of great magnitude. 26 


However, for more distant localities, where the runoff cycles coincidentally - 
may be radically | different, further clarification and the a of methods 


ofa adjustment are needed. 


This: paper _ offers suggestions for many interesting exploratory studies, 
especially for the younger engineers training in these fields who are afraid 


of “undertaking many laborious calculations. _ 


ALDEN Foster,” M. ASCE.- —The authors have a new 

a interesting method for classifying the daily-flow duration curves of rst _ 
a streams by a single ‘ ‘variability index,” together with the mean flow of the 
stream. At first glance, the variability index might seem to be related in some ~ 

way to the classification of skew probability curves by the coefficient of varia- 
- tion and the coefficient of skew, : as originally proposed by the late Allen Hazen,'” 
M. ASCE, and further developed by the writer.’ 15 _ Howe ever, the latter classi- 
fication was intended to be used in the study of selected samples of the , original 


Chf., Power Studies Branch, Div. of Water Control TV A, Knoxville, Tenn. 
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12 ‘Storage to Be Provided in Impounding Reservoirs for Municipal Water Supply,’’ by Allen Hazen, 
Transactions, ASCE, Vol. LXXVII, December, 1914,p.15389,0 
‘Theoretical Curves and Their Application to Engineering Problems,” by H. Alden 
Foster, ibid., Vol. LXXXVII, 1924, p. 142. Be 
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FOSTER: ON STREAM FLOW 


of actually ‘used was limited, s that the calculations be 
without an excessive amount of work. To compute the coefficients of variation 7 
_ and skew fora daily- flow record of several ye years would be entirely impracticable. - 
- Certain points related to the proposed method may need further clarifica- 


. a) The calculation of the variability ty index i is | independent | of the unit of 


° runoff used, whether cubic feet pe per - second or cubic feet per per second | per square 
mile, or flow i in terms of mean flow —because any change in in the unit —_ not 


tion: 


to shift the entire curve parallel to —_ til 


> 


(2) The variability index is a function only of the slope of the straight-line - 

“duration curve when plotted on logarithmic probability ‘paper. “Hence, the 
value of the index i As not changed if the curve is shifted parallel to itself on the 


(8) WwW Vhen applying the method fo: for a study i in n which t the ¢ data for plotting a 
x duration curve are incomplete, the following procedure should 


(a) ‘Assume a value for the variability index—I; 


(b) Determine a ‘value, Qi, such that log Q; = I (Q; renee a daily flow value, ue, 


aa expressed i in terms of the median stream flow); 
* (c) Assume a value for the median stream flow (the flow exceeded 50% 
the time), which m may be estimated from the mean flow, as: explained 
_ (d) Compute the discharge « exceeded 15. 87% of the t time as the product of — 
Q; and Qso (derivation of this percentage value being given subsequently); and 
— (e) Then plot the duration curve as a straight line on the logarithmic — 


per, ‘passing through the stream flow values for 15. 87% 7% and 


= 


Probability —To clarify the foregoing procedure, some explanation of 
the method of constructing the logarithmic probability paper is necessary. As_ 
- far as the writer - knows, 1 no adequate description of such a method se 
_ tion has been published. 7 Construction of this plotting paper ‘adem on the 


use of the “ “normal curve of error,’ ? which form: 


and is plo tted in in Fig. 4(a). The maximum ordinate of the curve is y. = o=h/Vr. 
The “probable | error” is r = 0.4769/h. _ The probable « error is also defined d by 
the fact that the area under the curve of error between the abscissas —r ‘and 
d +r is one half of the total area under the curve. The normal curve of error is 
known as a “frequency” curve. 7 If the frequency curve is integrated, the 


result is” a “duration curve,” or in the ‘present case the curve of the “prob- 
” as shown i in hen Plotted in in this the axis 
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FOSTER ‘ON STREAM FLOW 
| 
the ordinates becomes the ‘ ‘percentage- of-time”’ scale, ‘commonly in 
plotting duration curves. ‘As shown in Fig. the part of the 
curve between z and =-+ r is. represented 1 by the distance betw een 
25% and in ‘the percentage- time scale. 
4 
(a) NORMAL CURVE 
ERROR 
i 
| 
4 
q 
PROBABILITY 
INTEGRAL CURVE; 
| Fie. 4 4.—RELATION OF PROBABILITY ‘TO PRoBABILITY PLorrina 
if the abscissas, a, of t the he frequency « curve are expressed in terms of 1 r, the 
_ area under the curve from —a/r to +a/r m may ‘be found from values of the 
integral, as tabulated in . mathematical textbooks. This 
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| ‘FOSTER ON STREAM FLOW 


‘Tables of the probability generally give values of. asa a function®® of 
Qa or of” a, r. If P, is the value of the probability integral for a “residual” of 
(area from —a/r to +a/r), then will be the percentage- of-time i interval 
on the duration curve between values of (Qso — «<< and (Qs0 + a). ‘The symbol — 
- Qso denotes the value of the median, which is also the mean when the frequency 


= 


of two stream flow values, Q; and Q., differ from the mea ater a, and a2, Gay 


their percentage-of. -of-time range on the du duration curve (shown a as P; and P2 in 
Fig. 4(b)), or their values of total frequency, 1 may be found from the probability 
_ table corresponding to the values of a/r and a2/r. _ This fact is illustrated a 
Fig. 4(6) in which a,/r is assumed as 1.0; the corresponding value of = 
= 2.0, P2 = 0.8226, ete 
a - ‘The duration curve of Fig. 4(b) is plotted as a full line at natural scale. It : 
ean be reduced to a Pe line by a suitable change in the spacing of ll 


_— frequency taken from the probability table | is 0.5000 (= P;). Similarly, with 


paper’ 


equal to ai/d2, to 


let 


” (asi in the probability table), the duration curve | 


will plot i ina line. 
“4 Percentage X=0.4796 a/r Distance from 


TABLE 3.—Comuration OF PROBABILITY ScaLe 
of time | § | Ge.) 


“probability 


~ 


(hamming Distance from 0.1% to 99. 9% As 7.62 In.; from 5 50% 7 
Either 0.1% or 99.9% As3.81In.) 


lay 


x 


oo 


8 “A Short Table of Integrals,”” by Benjamin Osgood Peirce, Ginn & Co., New York, N. Y., 1910. - 


“*American Civil Engineers’ Handbook,’”’ Thaddeus Merriman, Editor, John & Sons, Inc., 
York, N. Y., 5th Ed., 1930, p. 79. 


4 


_ — For example, assume e that the distance on the plotting paper between 0.1% © ; 


and 99. 9% is to be d inches, and it t is desired to. determine the proper scale — 
distance betwe een 25% and 75%: P 


0.1791 
0.3708 
(0.4769 
0.5951 
0.9062 
0.9538 
1.1631 
1.4307 
(1.4522 
1.6449 
1.8213 


= 75% — 25%, or 0.50 
‘he 


Inc, 


X = ah = 0.4769 (6) | 
a? 
al 
= 
‘ 
8.87] 91.13 | 0.8226 2.0000 
5.0 | 950 | 0:90 
0.5 | 99.5 099 38190 | 3 
| >) = 


1.00: 1.00; and Pos = 99. 1.9% 0. 1% = 0.998, for which a/r = 4. 5817. 
seale distance from 25% to 75% will then equal (1.00/4.5817)d in. -Computa-— 
. 7 tion of other plotting positions on the probability paper is shown in Table 3. 
The values of X, corresponding to P, are taken from a published table’ 
f 7 * Percentage « of Time for Plotting Variability Index.— —W hen a straight line is 
Plotted on arithmetical (arithmetical scale for 
4 
E1.0 


Percentage of Times 


the at any any given of time 

that for 50% of time will be proportional to the probability integral, Pa, core 

responding to the assumed. percentage. Similarly, if a logarithmic scale is 

used for the ordinates, the corresponding difference between the logarithms — 


will b be proportional to The problem is to determine the percentage of 
time me at which the logarithmic ordinate of the : _assumed straight line has” a _ 


= 


| 
| 
| 
4 
_ Fig. Curves: Log Propasinity ScALE 


FOSTER ON STREAM FLOW 


95%. and computing the index as by the 
_ bata a more accurate res result will be obtained by assuming logarithmic — 


rain: 
‘proportional to the corresponding values of the probability integral, computing 
corresponding values of the index, and then determining the 
of time or probability (P.) that would have a value of a/r equal to the index. — 3 


ara index > 


Locus of means | 


Stream Flow in Terms of Median Flow | 
o 


an 

q 
= SS 
00 a 
_ The computation is st shown in Table 4, in which Col. 1 denotes the range in | 
"percentage of time (= 50% : 4 P,/2), and the data i in Col. 3 are selected from 7 ; 


available probability tables. 28 The index i computed with 


in Col. 5, thus: The in index from from Table Table 4 
14660. 


Percentage of time 


15 to 85 70 7329 2. 3617 

0.30 

0.0847 


| 
A 
| 


> 
> 


= 


The corresponding value of X = 1.4660 times 0.4769 = 0.6991; and a 


‘id corresponding to X = 0. .6991, is 0.67718. Then, the percentage of time, for 


the computed index, is +50 = 16. 16.14%. The discrepancy, 


5.—CuEck CoMPUTATION (OF VARIABILITY 


0.8228 0.67700 
0. a | (0.26884 


= 


0.3365 0.11323 
0.1931 


0.03728 
0. 


organ 


20150 


this value and 15. 87% as as given by the authors may be due to the omission of : 
‘the decimals, 
series of duration curves with va h values 0 of the index from 0. 1 
1.0 is shown in Fig. 5. For con- 
venience in plotting, the 50% 
median stream flow value is 
sumed as 1.0 ‘(logarithm = 0.0) 
For an index of 1.0, the line will 
pass through a value of 10.0 
(logarithm : = 1.9) at 15. 87% of 7 08 
‘The other lines will | = 


arithmetical ordinates, at the same 


percentage of time, equal to their — 
respective index values. 
show the true relative shape 
of these duration curves, they have 
been replotted in Fig. 6 on arith- 
metical scales. A check of the 
accuracy this: ‘method is’ given 
“in Table 5, in which the values of 
flow (Q) were scaled from the line 
plotted for the index = 0. 50. on 


the case of Table 4, the vari- 


‘ 
percentage of time—— 
Ratio io of Median to ean an Value 02 
Stream Flow.—Unless the dura- 


4 curve of "stream f flow has been 0.1 


= 


Ratio 


oe 04 06 O08 10 


actually calculated, there is no 


although the mean stream flow is generally known or may be estimated. — The 


method of constructing the duration ¢ curve ‘Proposed the authors requires 
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OSPINA AND STREAM F FLOW 


“a the u use of the median, and consequently it is 3 necessary ry to establish a relation 


tempted to ‘establish by integration of the. plotted 
curves. a The results are shown in Fig. 7 which also indicates the > percentage - 
_ of time at which the duration curve crosses the mean ordinate. erie 
S. Osprva® AND G. Jun. ASCE. —The “variability index” offers 
an excellent method for the hydrologist to compare conditions in different water-_ 
= = _ Its computation is simple and it does quite accurately represent the — ; 
greater part of the flow variations. To facilitate plotting adjusted curves 
- fined by the mean flow and the variability index, the writers suggest the con- 
struction of a graph which shows the of time during which the mean 


flow is exceeded, as a function of the index For indexes 


ae 


ae 


OVS 


oO 


FLOW RATIO (Qm. =1.0) 


— 


> 


7 


2 


VARIABILITY INDEXES 


(0.10 
0.00 
10 
20 20 30 40 50 60 70 90 "100 


8.—ApsusteD Unit DurRATION CURVES; Fic. 9.—Apsustep Duration CuRVES FoR 
PercenTace oF Time Durina: —— Durvenane | VARIABILITY INDEXES 


“greater than 0.35 ‘only | a range can be indicated, with sufficient p precision fo 
_ practical purposes, because the average is affected by high discharges occurring 


during a smé all percentage of the time. Fig. 9 shows the adjusted duration 


curves used, which are straight lines on eaten mic Probability p paper. wees 


4 


30 Parteer, Olarte, Ospina, end Pay4n, Engrs., Cali, Colombia. 
31 Associate, Olarte, Ospina, Arias and om. Colombia. 
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“4 OSPINA AND TAMA ON STREAM FLOW 

| ‘The pub published data reveal no definite relation between mean unit flow 

(cubic feet per square mile) and | the variability index, except for a certain 

- similarity in rivers situated in the s same e part of oe United States. It is noted - 


than 1 000 sq aan on and that the authors’ conclusions in regard to the — 
a ' bility iadiex incline to refer to small rivers. Both common sense and study of _ 
the tables confirm a definite i in 
"fluence of the size of the TABLE 6. 
shed on the index. _ The values a 
in Table 6 have been deduced 4 
the data presented by 


Drainage area No. oat Most Maxi- | Mini- 


(square miles) | stations frequent mum mum 


0.1 to 1,000..... 0.40 t00.50 | 1.17 | 0.14 
t ‘1 
authors. It is apparent that 87 | 00.38 | O98 | 


extremely “high “indexes occur 


a that there is a general tendency toward lower values as the drainage area 


the Cauca basin in W estern Colombia the following information can 
bes given: Flow variations on the main river and most of the tributaries follow 
the same trend, regardless of the size of the watershed. An index of 0.23 is - 


: typical. The mean annual runoff is ‘exceeded during 40% « of the time and it 


. averages 1.8 cu ft per sec per sq mile. In contrast to the uniformity of the 
- variability index the unit runoff varies considerably even in watersheds with - 
; similar rainfall characteristics. 3. The total drainage area at the lower end of 


the so-called Cauca Valley is 9,000 sq miles, the greater part of which is located : 


in rugged country. — Soils vary from impervious silts to pervious alluvial de-— - 

posits in the flat regions; in most of the mountainous regions the soils are im- __ 
pervious with ‘grass, bush, and forest cover. There are two dry seasons and 
- two rainy seasons s and the annual rainfall is about 40 in. ares 
After years of work in irrigation and power planning, | the writers were =. 

- surprised to find a variability index so low in comparison to that for streamsin 
the United tates 
The authors recommend their procedure for studies where the available 
toca records are meager a and they rightly state some of the factors swhich — 
7 limit its use. © _ The annual mean flow is to be determined from short-term flow | 
data, from comparison with other hydrological records and drainage basins, — 
from temperature records, ete. _ The selection of the variability index is to be 
governed by geological conditions and surface storage only, whereas rainfall — 
distribution and intensity, vegetation | cover, and arrangement of tributaries 


— 


_ are acknowledged to be of influence to an unknown degree a oe 
7 Actually, the number of factors is such that no accurate predictions can 
be made unless flow records exist ‘covering a certain p period. ¥ To the writers’ = 
knowledge, even short-term discharge data are likely 1 to be more useful than 7 
coefficients derived from conditions in other regions. Attention As called to 
the fact that geologic conditions may be defined for restricted are areas, yet i in large _s 
basins the problem _ will become rather complex. Preliminary estimates 


‘stream, flow i in countries are > hampered by deficient hydrological 7 
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> duration curves, are on mean flow and | have the 
_ being dimensionless, as well as chronological records from other rivers, are 


records should be checked, since systematic errors in high discharges, due to 


~ deficient section rating, will affect the variability index. | 


=| 


= 


FLOW 
© 15 YEAR RECORD 


©) YEAR RECORD = 


@ YEAR R RECORD 


© 6 MONTH RECORD 


30 40 50 60 70 80 90 
10. oF Duration Curves ror DiFFERENT RIODS ON THE Cauca RIVER 


experience, since both mean flow and variability index change w 
the period considered. 
‘ river (Cauca) for different periods, related to ‘the 1 15-year average of Qu = . 
‘The: mean flow of the other periods is 3 indicated by circles. Iti is evident that 
* “fndex” method of adjusting duration curves to straight lines on loga- 
rithmic Probability paper should * records. It 


‘OSPINA “AND “TAMA ON — 


The first guess will teflect Laer 


- valuable to deter mine flow variations: W ariability index), annual mean dis- 
_ charge, , and duration of droughts c on the river studied. Mass curves of rainfall — 
os will be found useful to analyze the general trend. . The accuracy of stream flow : 


adjustment ‘short- term. records | to fit : average conditions requires 
vith the length _ 
Fig. 10 shows the actual duration curves of the — 


— 


| 

ig 

fo 


¥ 
the drainage basin. 
estimating method i 
established gaging stations win Provide more Tenable Ulder 
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KAZMANN ON ON FLOW ] 
- ‘is necessary to examine thoroughly a all pertinent relations until a a representative 
duration curve can be plotted. 
a Finally, it should be remembered that: most hydraulic s studies will require a an _ 
estimate of critical conditions during a period of severe drought. These are 
not fully defined by a duration curve representing a long-term record, and they — 
have to be determined | by a design hydrograph based on actual « Cachan data 


and adjusted to the most critical conditions. 

‘The writers’ experience on the hydrological studies in the Cauca Valley 
4 


. _ shows tl that, when only short- term records are available, the most difficult task 
| the determination of the minimum flow hydrograph—w hich represents 
: minimum flow itself but also so the hydrograph for the entire ¢ critical dry season. 7 
; a This hydrograph i is of the | greatest importance in planning irrigation and pow ower 
“roe where good storage reservoir sites are not available. _ ‘The ratio. of 
minimum flow to average flow for a period of from 3 to 15 years varies from 

| :2.3 to 1:5.8 for the tributaries of the Cauca River, all with practically equal 
variability indexes. The Cauca itself has a ratio of 1: 1:3.6 6. On the other hand, 
a the characteristics of the dry : season hy ydrograph vary’ vary very widely from one 


| this: indicates the complexities and uncertainties 2s of hydrological extra-— 

polations. — _ The writers consider that the greatest usefulness of the variability — 

index, w here few records are available, is to help develop t the judgment of the 

7 engineer -and thus aid him to make a better. guess of the 1 mean and the minimum — 

_, The principal value of the paper is considered to be t the development. of a 

useful statistical investigating procedure and the analysis of a considerable oo 

~ volume of flow data from North American rivers. It is then hoped that other 

will contribute information from different regions and nations in 


order to further the study. of the factors which affect the index, 

RAPHAEL G. KazMANN, Assoc. M. ASCE.—The authors w ‘ould do a a great 
service if they would enlarge on the use of their variability index. For example, 
- the method depends on the proper ‘selection o of the v variability index as outlined 

| 


on 


= 


in Table 2. This table translates qualitative descriptions of geology, hydrol- 
= and topography into quantitative terms. _ Apparently a a rather aor with. 


WwW hat permeabilities and what niennane of the drainage area are meant 
; by the various subheadings 1 under “bedrock and soil?” What slopes and 


percentages of drainage area ai are ‘referred to under “relief?” Likewise what, a 
quantitatively, is the average permeability of moraines? Under “lakes and 
swamps” what percentages of drainage area correspond to ‘ “moder 


: T hen, there is the problem of determining the mean annual runoff of the — 
stream. ‘Using the method given by Don. Johnsto one an and 


‘Engr., , Ranney ‘Method Water ‘Supplies, Columbus, Ohio. 


> 
i 
1 
2 
— 
— 
| 4 
! 
| > 
if 
= @ 
— 
: 


IN AND HUDSON STREAM FLOW 


illiam Perry Cross, M. ASCE, t the of stream flow re record needed 


estimate the mean annual flow to various degrees of accuracy, for three Ohio 
streams, is shown in Table 7. These streams ns have been n gaged continuously 
| TABLE 7.—LENcTH or REcorD NEEDED” | 
Licking...| Toboso, Ohi 672 10 
Springfield, Ohio | 48 485 53 = =. 10 | 
1923, Ons an ungaged stream, how many years of gaging do the authors 
st, er what degree of accuracy is desired, i in determining the mean annual a 3 s 
o Bi. The writer would like to withhold judgment on on the proposed ‘method of ; b 
developing flow duration curves until studies are made on watersheds whose | 
physical characteristics are known and where long- -term gaging records are si 
The magnitude of deviation between the computed duration. curves . ¢ n 
a and the actual duration curves would show conclusively | how applicable Table 2 _ a Z 
to field conditions. Naturally, all the suggested studies would have to be 
made w without reference to actual long-term gaging records if valid comparisons th 
are to to result. | 
Brurn™ anp H. E. Hupson, M. ASCE.—The | proposed. method 
=; of evaluating stream flow variability is useful in applying data 
streams with long periods of record to comparable streams for which 
There have been a number of developments: of small Illinois L 
i: for municipal supply | purposes. Iti is expected that there will be an increase i - R 
in the use of smaller watersheds for impounded supplies to minimize sediment _ p 
damage and spillway costs. . Periods of zero flow exceeding 5% of the time a ia th 
are often encountered in such streams, making the method of this pa paper inap- 7 7 
plicable. Similar situations are found in flowing on highly permeable 
valley fills and in arid regions where rainfall is infrequent. The method i is: 
worthy of modification to make it applicable to intermittent streams. 
“a An attempt was made to compare t the stream flow in Money Creek ae | a: 
Bloomington with other gage records in central Illinois. For this purpose flc 
stream flow variability indexes were computed. An immediate difficulty made 
‘itself evident i in that Money Creek had a zero flow 547% of the 12- -year q bj 
"Period of record. The variability index became infinity. 
An inspection was made of the stream flow records from Illinois gaging _ : 
; "stations, and it was found that this problem existed for five watersheds in the : ' 
‘Elements of Applied Hydrology,” by Don nd William Cross, Ronald Press Co., 
Eng. Asst., Eng. State Water Survey Div., Devt of Registration Education, Urbans, 
Head, Eng. Sub-Div., State Water Savey Div., Dept. of Registration an and Eduestion, Urbana, Ill. 


+ 
: 


group from 45 sq miles to 500 - miles. ' The largest of the gaged areas sid 
zero flow more than 5% of f the t time is Mazon . River, which drains a watershed 


‘aad exists i in the two watersheds (areas in square miles): 


Slowly to very slowly. 


through cooperation — with the 


United St States Geological Survey, 
are shown i in Fig. 11. 


an In view of the inapplieability 

of the authors ’ method to 

streams in which there is zer ro 


of the average subsoil permea: a- 
bility. in Illinois, whereas the e 
Mazon River subsoils” are con n- 
Zero. flo ow was reported for 
Mason River during 7. 06% of 
the ‘period ‘of record, which is” TTT 
9 years, s. this period the run- 
from streams in the same 
region vas slightly above normal. 
Lake Bloomington and Mazon £ | 
2 


as 


flow more than 5% of the time, . 
it is suggested that the varia-— 


4 
Calculated values for the streams Fig. 11.—Drscuarce Duration Curve, Money Crezk 
AND Mazon Rt RIVER IN ILLINOIS 
discussed herein, using. the per- 


- centage values and procedure proposed by the authors are as follows: ———— 


‘Physiegraphic Divisions of Illinois,”’ M. George E. Ekblaw, and Leland Horberg, 
The of Geology, January, 1948, p. 16. 


87 *Principal Soil Association Areas of Illinois,” ‘Dept. of Agronomy, Univ. of Illinois Agri. Experiment 7 7 


‘| 
BRUIN AND HUDSON ON STREAM FLOW 
— 
Mon eek drains 9 watershed of 45 sq mi which lies between two 
| oil. 
] 
— 


—_ 


Percentage of time of zero flow 7 18.540 


_ This index would be accompanied by a dates of = proportion « of time 


4 during which there is zero flow. eee 


ee EW. Luu.*—A particularly interesting feature of this ingenious index of 
jet 7 variability, not pointed out by the authors, is its s comparability in in Meaning toa 
—_— a coefficient of variation—that is, it represents n not variation, as such, but variation — 

in respect to the mean. Since a constant logarithmic deviation corresponds to 
a constant percentage change, computing the standard deviation of the logar- nue 
7 1 ithms of the discharge gives an index of deviations in m percentage of the mean, — 


— 


as is done directly by the coefficient of variation, 
r 


0 


An undesirable feature of the index is that it cannot be 
mean flow. However, because of the normality of the distribution of en 


— in throughout most of the range, the index when used 


_ utility of the index would be increased if estimates of discharge at the 1% — 

— 99% frequencies were appended to it. _ These values, by indicating the nature — : 
of the curvilinearity of the extremes, would be of importance to those interested =! 
either in flood or in minimum flows. 
An interesting 1g point is s the relation of the variability indexes to the drainage 
area. The authors state that streams from large watersheds will tend to have 
_ indexes than streams from small watersheds because of less synchroni- 
gation of flow and greater channel storage. - Surprisingly, a a tabulation of the ~ 
area and 1 the variation of 22 streams, | listed by t the 2 authors (each of ¥ of which were 


gaged at 2 or more : stations), reveals that this is not the case (see Table 8). 
_ Seemingly, the indexes should decrease in the downstream direction. _ Accord- 


U 


ing to the tabulation, this occurred on 11 streams, whereas on 8 streams the 
indexes became larger, ‘and on 3 streams (Table 8(a)) they remained constant. 
Closely associated with a decrease i 2 in the variability indexes is an increase in 
— the mean flow; and with an increase in the indexes, there is a decrease in the 
; mean flow (see Cols. 4, Table 8). _ Since the indexes reflect variation in respect 
to the mean, these relationships permit the e speculation. that either variation 


remains e constant, or, 1 more likely, that it changes i in the same direction but less 
= ‘rapidly than does mean flow 
After considering several factors which affect s stream flow variability, 
=e authors concluded that geology and soils, and surface storage had the major — 
Although there is some evidence to support: this conclusion, it ‘is 
hardly sufficient to warrant the finality the authors give to it. . In a study | of 
this kind, it is extremely doubtful if _— having an lenge bearing “a 
stream flow w variability can be isolated. 


One reason for this doubt is the posed the w range 


areas. On small drainages, vegetation and weather may be important 
a Div. of Forest Influences, Forest Service, U. S. Dept. of Agricul nm aati D. C. 
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LULL ON STREAM FLOW 


larger areas. Geology and soils could peared most factor on on 


areas possessing a variety of ‘different types of vegetation. Similarly, 


TABLE 8.—RELATION OF ‘STREAM FLow VARIATION TO Size or DRAINAGE 


Drainage ndex¢ Mean | Sta- | Drainage Index* Mean | Sta- Drainage Mean 
4 area® flow? | tions | area’ flow? | tion? flowé 
(2) @ @ (2). (2) ) | (4) | 
a) VARIATION REMAINS © ARIATION DECREASES WITH | ees lilies 
4 WiTH INCREASE Con INAREA (c) ‘Vastarion Incapasee 
Muskingum River, Ohio Wapsipinicon River, lowa St. Croix River, Wisconsin 


ya 550 | 0.19 | 0.761 
2'390 0.24 | 0.705 


36 | 5,982 | 0.48 | 1.021 | 29 | 1,060 | 0.64 | 0.442 | 
7411 | 0.48 | 0.791 2'300 0.564 
Lamoille River, Vermont on 


51120 | 031 | 0.583 
(Table Sandusky River, Ohio (Table 1(d)) | 5,930 0.30 | 0.566 
1(9)) 812 | 0.20 5 
7 it 43 River, Ohio (Table ua) | 29 | 2240 | 038 
wits > O88 | | 2.340 | 0.32 | 0.762 
IncREASE mm Anza Miami River, Ohio (Table 1(d)) 4, 870 | 0.42 | 0.640 
q Ss - 
ea, 39 3639 0.46 | 0.966 | ¥ 550 | 0.76 | 0.763 
| 5,490 | 0.64 | 0.278 | | | 7 
6,818 | 0.63 | 0.267 | Ashuelot t River, New 
6 13,200 Maumee River, Oni (Table 1(d)) 
| 13;900 | 0:50 | 0.382 | © E 56 | 0.744 
Skunk River, Iowa 420 | 0:45 2 0.79 | 0.671 
(Table 1(a)) | 420 5 
| 320 | 0.94 | 0.350 Conne ecticut River, — River, Ohio (Table 1( 
10 | 2,800 | 0.53 | 0.495 “(tab |, 333 | 0.56 | 0.883 
(Table 8,000 | 9. 29 | 0.73 | 0.606 
iver, lowa (Table 1(a)) 4 
4 13° | 3,230 | 0.49 0.457 Delaware York and | 23 1,398 0. 000 
12/480 0.38 | 0.489 2 "439 0. 47 1.007 
Cedar River, Iowa (Table 1(a)) | New Jersey (Table 1(/)) 4 
‘15 | | 845 | 0.39 | 0.341 | 4 7 | 875 | 0.47 | 2.331 Mad River, Ohio (Table 1(4)) 
ig | 1,660 0.38 0.407 13070 | 0.44 | | 485 | 0.29 | 1. 053 
| 6,640 | 0.36 | 0.442 | 6,800 | 0.35 | 1.690 | 29 | 632 | 0.30 | 1.050_ 


@ Cols. 1, 2, and 3, Table 1. * Col. 4, Table 1. © Col.6, Table 1. 4 Col. 7, 


an entirely different character than larger drainages including a a diversity. of topo- 
graphical features. Thus, it is difficult to select factors of | major or importance e 
watersheds of different sizes are considered. = 
_ Furthermore, unless each of the several factors is given equal ‘consideration, — 
is doubtful if any conclusion is justified as as to relative importance. At 
the present time, lack of data on all factors would prohibit such an analysis; : 
but. it i may | be that if the authors had given more attention to factors considered - 
~ of minor influence, a different conclusion would have been reached as to relative 
For instance, from the authors’ conclusions the iaiiatiaiin can be made that 


_ vegetation is not an influential factor. A comparison of the variation of streams 7 


fac or droughts, or density 

— 

| 

| 

| 
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from forested and nonforested areas 2as of size in Iowa and Wi isconsin does 
not support this view (see | Table 9). _ Furthermore, it should be emphasized 


that the relatively s small ‘variation fr from areas of greatest t relief, in I Iowa and 
: North Carolina, is perhaps due more to the fact that they are mainly forest 7 ee 
; - covered than to the reasons given by the authors. If weather were truly con- _ 


sidered a factor, perhaps the authors should have the effect of the 


TABLE 9. —ComParison OF VARIABILITY INDEXES oF ForRESTED AND 


NoNFORESTED ep AREAS (Sq Mixes)» 


ndex 


— 


0.94 

224 Lime Creek... 0.63 


Sugar 

‘Peshtigo. . 


“Mean variability (Table (a)) .........| 0.39 | Mean variability (Table (a)) ......... 0.49 


states. In | addition, the fa factors deemed in important by the authors, when con- 

sidered singly and in respect to the area involved, are not of major importance. : 
Surface storage can obviously be an important factor only in the “4 


: ; limited areas in which it exists. — Geology and soils were found to be important 
in n certain ; areas, but not in enough a areas to justify selection as a ‘major factor. 
_ Perhaps a more reasonable conclusion n would be th that the > relative ir importance E 
ee of the factors cannot be determined for as large | an area as the one concerned. 
All the factors may be considered as having a bearing on variability; each is 
: present t to a degree in all drainages ; in certain n drainages one or two factors may 
be of particular importance because ‘they are present in sufficient intensity to 
mask the effect of the other factors. However, because of the variability of the a 
~ factors over the particular area, in addition to the variation in size of the drain- 


tends to intensity the effect of the ‘more variable factors on the 


smaller areas, and the fact that evidence is not available to permit a thorough 
= - analysis of the effect of each factor, the data that the authors present do not — 


allow an evaluation of the relative importa nee 0 f the factors affecting stream — 
flow variability, 


-* 

a 

— 

Pine | 0.32 | Milwaukee 661 | 0.48 
7 
| 
fl 


7 


D. MrrcHett, _ASCE- —The authors } have made a sound 
i 7 and comprehensive presentation of stream flow duration. Their - use of log- 
-arithmic probability scales is particularly commendable. 
7 J Instances are rare indeed i in Ww hich the logarithmic probability plot of the 
4 flow duration curve will yield a straight line throughout the entire range of 
flow. Forsuch rare cases, the ¢ curve is adequately described by the two param-— 
eters by the -authors—namely, the variability index and the mean 
discharge. — As they have pointed out, the position of the curve must be > estab- — 
lished by trial, but that. task is not particularly difficult. On the other hand, 
in all cases, with which the writer is familiar, there is substantial curvature — 
: near one or both ends of the logarithmic probability plot so that the two 
- parameters : are not adequate to define the curve. “- This fact has been recognized 4 
_ by the authors, as evidenced by their statement (under the heading, “Use of 
the ‘Variability Index”): “ ‘The ends of this line for a short distance on each 
end will have to be altered to match the best possible estimates of maximum = 
- A project on which the writer is engaged involves the flow duration curves — 
7 q for 15 of the Illinois stations used by the authors. The basic data differ in in = 
one major respect—namely, the records have been extended to 
ki complete. water years through 1945. One of the stations, Big Muddy 
River at Plumfield, has been selected to. show the between the 


from the station data. 


7 and the mean discharge per square mile (q) were known to as. given 
Table 1, Cols. 6 and 7, 1.01 and 1. 021, respectively. From published 
® data,” it was determined that the extremes of discharge prior to September 


- 30, 1933, were 16,300 cu ft per sec and 0 cu ft per sec (the equivalent of 21.6 = 


4 cu ft per sec per sq q mile and 0 cu ft per sec per sq mile) and it was inferred — § 
that the time of no flow was about 9 days in 19 years, or a little more than 0. 1%. 
The resulting synthetic curve appears as the continuous line in Fig. | ; 


he In determining the slope of this line, it was necessary to ‘modify the statement 


of the authors (under the heading, “Use of the Index’’): 


“The s shape of the duration curve can be obtained by drawing a straight- 

3% line duration curve on logarithmic probability paper with a slope such — 
that the ratio of the dischar; ge exceeded 15.87% of the time to the Soon > 
a exceeded 50% of the time was equal to the variability index selected.” ai 
Since the selected index is 1.01, it appears | that a line of such slope would be 


nearly horizontal, and not at all representative « of the flow duration curve . 


Obviously, the ratio of the discharges should be equal to the. antilogarithm 
; of the selected index, or, in this case, 10.23. The line was drawn on this slope : 


7 “Mor on the first trial, Ws vas made to pass through the point | given by J. H. : 


Morgan"—q = 1,01 and percentage of time equals 25.8. The data on extreme 


“Surface Water Supply of the United States, 1933, Water- Supply P Paper No. 745, U.S.G.S., ., Wash- 


4 ‘‘Flow-Duration Characteristics of Illinois Streams,” by J. H. Morgan, Transactions, Am. Geo- 
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- discharges offered no reason for curvature in the lower end of the line, but the 
eo end obviously must become asymptotic to the horizontal q = 21.6, and — 


charge, to the use se of. addi- 
tional record, is largely responsible 

; for the downward | shift from the | 


Bia. 12. —Dunasson Curves ror Fiows sor synthetic to the actual ct curve. 
Bra Muppy River ar = a percentage basis this shift is a 


maximum at about the 70% dura- 
tion point, where the synthetic curve a disc charge ‘of more than 

7 double that of the actual curve. ~ Had this particular point been the ne subject 
of investigation, the use of the _ synthetic curve could not be recommended. _ 4 
_ At other points, particularly in the region of the mean discharge, the agreement 


between the two curves is very satisfactory. It should be noted that the 
- assumed shape of the synthetic curve has been helpful i in obtaining the good 


the line was drawn accordingly. ¥ Step integration of the area under this line a 
“yielded a mean discharge of 1. .25 cu cu ft t per sec per sq mile, greatly i in excess of : 
‘per sec per sq_ mile, which ‘indi- 
cated that the curve must be 
=: moved leftward to the position as 
‘For comparison, the duration ‘ 
curve as developed by the writer, 
z SS based on ac actual r records for the 35 
ater years, from 1909 to 1912 and 
35° +H from 1915 to 1945, is shown as the 
dashed line in Fig. 12. As indi- 
‘ated in Table 10, the variability 
‘index for the e longer period i is 1.08 
the mean discharge is 0. 891 
5 £2 eu ft per sec per sq mile. The 
— duration of the mean discharge is 
23.5% of the total time. 
The” decrease in mean dis- 


ime in Percentage of Total Period 


agreement in | the Tegion near the mean discharge. Since the parameters fo: for = 
the synthetic ¢ curve are on the | same data as used by Mr. Morgan, 


inaccuracy by shifting the synthetic curve to pass through the point g = 1.01 7 = 
and the percentage | of time equals 21.8. It is believed that the improvement _ 
7 _ thus obt obtained i is only coincidental; in other cases the shift necessary | to nica coll uy 
sate for inaccurate shape may be ‘such as to. move the synthetic « curve away 

from the actual curve in the region of the mean discharge. 

_ Furthermore, the curvature in the upper part of Fig. 12 is opposite to that 


mentioned by the authors: (unde the heading, ‘ “Use of the Index”). 
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This is due to the fact that Big Muddy —— is — by extensive food 


_ plains which, by their storage e effects, materially reduce the magnitude | of the 


. _ The foregoing e: example t benefits, of course, from the fa fact that very pertinent — 


per were assumed for the index and the mean un discharge. Had it been ne neces- 


TABLE 10.—CHancE IN VARIABILITY INDEX AND Mean Discanos 
SELECTED ‘SITEs, To AppiTion To LencTH or Reconp; 


* IN (SEE | TABLE WW) : 


>= 


Lane-Lei Mitchell Ratio | Mitchell 


Pecatonica 
Algonquin 
Lyndon 0.30 
Momence 
Seville 


— 


La Moine . Ripley 

Sangamon Riverton 

Sangamon Monticello | 
| Vermilion Danville _ 

Macoupin Creek.| Kane 

South Fork Kincaid 

Kaskaskia Vandalia 

Skillet Fork ....| Wayne City 

Big Muddy Plumfield 


| 


an 


oornd 


Forman 


| 
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7 sary to estimate these values from data given in Table 2 or from records of | 

* ‘hear- ‘by watersheds, any less appropriate e estimates would have been 

>. The writer has had no opportunity to test the data given in Table 2. 

Some information may be provided, however, on the variation of the index 
at a given point ona given stream, as a result of 7 varia 1 the period of — 
record. As stated, new duration curves have been prepared for 15 of the 
described by the e authors. These new curves include data for all 
complete water years prior to 1946. i values “of the mean discharge and 
of the variability index, computed by the authors’ method, are shown in | Table : 

stations, the index value has. changed appreciably, and, in 

_ general, the index value has increased while the mean discharge has decreased. 
These fifteen cases are too few to warrant firm conclusions as to trends; but 
they do raise a question as to whether an increase in the value of the ae 
index may not be expected with an increase in the length of record - “4 _ 


OwnseEy,® Jun. ASCE. —This paper should help to stimulate the 
“interest of hydraulic engineers in the possibilities of a powerful tool for the 
evaluation of hydrologic data, which heretofore has not received the attention nit 

deserves. With some notable exceptions, hydraulic engineers have not kept’ 
- ‘pace wi with, or taken proper advantage of of, the developments i in the mathematical 


 seience of statistics. ai authors have introduced, to the analys sis of stream 


43 Asst. Prof. of Civ. Eng., Univ. of Tennessee, Knoxville, Tenn. 
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flow data, one of the characterizing numbers used by statisticians as a measure — a 


_ of the nature of the distribution of data. 


introducing their variability index Sie the heading, ‘ “Previous 
‘Studies of Duration Curve Properties”) they s state that: 
“Many civil engineers are familiar with this from its application in triangula- _ 
_ tion and base line measurement, where the standard deviation is known as 
the ‘root mean square error.” 
2 They might also have stated another analogy that is perhaps even more familiar +. 


civil engineers—that of of the mathematical and physical similarity to the 


moment of inertia and the radius of gyration yn of areas. - Inc characterizing s stream 2 
flow data, the engineer is interested in the spread, or dispersion, of values, — 
- rather than in the probable error; and this is essentially what the radius of 


gyration measures for the cross-sectional area of a column « or other structural i - .o 


Without having seen the original, n more detailed paper paper filed in the Engi- 7 
neering Societies Library,®* the w riter believes that this paper has suffered in 
_ condensation from the omission of background | material. . The explanation of 
the method of constructing a synthetic frequency curve from the variability 
_ index and the mean flow, in particular, is not cl clear, we iii 7 
Iti is wondered why the authors chose to depart from the accepted definition — 


a “standard deviation” in defining their variability index. ‘defining 


in which -X is an . individual variate; 2g is the mean value; and N is the total 
number of variates in the population | group, w which, in the case of daily stream 
-s flows, would mean the number of days from minus to > plus eternity. W hen the : 
4 a. calculation of standard deviation is based on a “sample’ ’ consisting of n days — 
7 of record, a fair approximation is obtained by applying a correction factor 
WVn/n — 1, making Eq. 7 equal to da 


The precision of the number so. obtained i increases with d the number « of f days, . 
——— used in the » calculation, and d it will be noted that the correction on factor ap- [ 
In effect, the authors have ap the correction in 8 by u using 
the divisor 9 for their calculations based on the 10 values picked from a fre- 
quency curve, which is an overcorrection 1 (as will be explained), — 
use of 10 numbers: picked from a graph i is questionable. lt seem 


better to group the raw data into equal class intervals and to use midvalues a 


these intervals, with corresponding frequencies, , for the calculation of the index 
‘These intervele should be chosen with regard for proper statistical ——— 
- tion of the data, but would surely be such as to give more than 10 groups in 


29 W. 39th St., New York 18, N. 
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case. At the sa the same the arithmetical simplifications made by” 
the use ‘of. -equally- spaced intervals would reduce the labor involved a 
minimum. Iti is doubtful that the calculation (suitably made) of the index 
4 from the raw data would require any more time than a computation by the 
method of picking 10 from the curve e as the 


=e 


again to the use of 9 as a divisor i in “variability 
index,” this would be the proper number to use if the index were based on only | 
10 values of raw data, except that, if the index had been based on only 10 days yo 
- of record, it would be statistically worthless. — The authors have stated that — 
the minimum period used in the compilation of their tables was 10 years, ' which © ? 
4 for daily flows means a minimum of 3,650 variates. — Would it not have been 
more nearly correct to divide by 10, the number so obtained bein; g theoretically 
to a correction factor equal to 3 1650/8, 649, Ww rhich is negligible? 
principal value of paper in these studies would seem to 


—the proof being that plotted on tin 

: paper yield a straight line. " Apparently the authors concluded that the actual 
~ values of stream flow do not follow the normal law, but that the logarithms ¢ of 
q these values do follow it—hence, the use of rather than actual 


stream flow in calculating the index, 
a In any problem for which it can be assumed that the variates (whether raw 
‘data or their logarithms) follow the normal | distribution curve, it is a 1 simple 
~ matter to compute » values for a frequency curve with percentage of time | as _ 
argument (called a flow-duration curve by hydraulic engineers). For this 
_ purpose | it is necessary to have only t the m mean | v ‘alue and the standard deviation | a 
4 (known or estimated from tables such as those presented i in this Lancig and a 
ever, one e characteristic of the normal Gitetvation | is ; equality of the mean and 4 


= the median. For stream flow it would be necessary to assume that the mean _ 7 


flow was the same as the 1 flow equaled or exceeded 1 50% « of the time, and this 


is not usually found to be. true. _ Perhaps the authors’ method of ' constructing b 
If so, it is 
that farther ape will be given in the closing discussion. 
‘The numerical value of the standard deviation as applied to ‘the flow of 
a: stream re reflects not only the variability of the flow but also its magnitude, 
The great value of date 1 such as presented i in 1 this paper should be in the possible - 
use of a stream as a “‘model,” from which to construct or check a flow-duration 7 - 


g 


curve for some comparable st stream for which data are scarce. | To eliminate the 
7 effect of size from such a ‘comparison, it seems that the ratio of the standard — 
deviation to the mean flow (the coefficient of variation) should be used, rather | 
than the standard ‘deviation itself. It should be e understood, of ponents that 
: 4 in order for 2 streams to be comparable, their sizes could not be greatly different _ - 
of the order of scale ratios commonly used for other | purposes 
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hn Cotton, M. 

ASCE, the fact that the determination of reservoir capacity needs. 
as much or more attention than most other current problems 
ban to which engineers habitually give their attention. The paper provides : addi- a 
tional techniques through which such capacity studies can be made. Nearly 
all countries possess rainfall information from which reasonable estimates of - 

mean annual runoff are possible; but knowledge of the variation of the daily — 
stream flow that determines reservoir capacity is ; much more limited. This 

= situation is particularly true in underdeveloped countries where knowledge 
of minimum s stream flow i is especially needed now because of current emphasis — a 
on_-wide-seale river basin planning. ‘The paper deals with | two. problems: 

Ta the development of an adequate index through» which daily stream flow 

variation can be expressed; and, second, the determination of the factors. from | 

a knowledge of which its quantitative value can be estimated. ee 
“ 7 The utility of arithmetic averages for summarizing a mass of data is com-_ 


“monplace. Hydrologists guide their judgment through their knowledge of 


mean runoff per square mile, mean rainfall, mean flood discharges, and so - 
forth. Knowledge of these mean values helps prevent mistakes when the 


engineer is dealing with ‘unreasonable observations or fragmentary data. 


——————— the other hand, acceptable data on the mean values of stream flow variation 

_ are few, perhaps because of the lack of a generally y accepted variation coefficient — - 

through which past studies could have been summarized. In placing before 

: the profession their logarithmic ‘‘flow variability index,” and in giving its 7 
ie - value for some 200 American streams as computed from hundreds of thousands © : 
of observations, the ar authors provide | hydrologists with type values for a 


= physical characteristics are w ell and hence provide the en; engineer 

—— Long ago Francis Gi Galton pointed out it that if, instead of characterizing t the aa 

eee distribution of sets of observations of natural phenomena such as : 

runoff, 1 rainfall, and so forth by standard deviations, they be be characterized by 

7 the standard deviation of the logarithm of the | observations, ‘a much better 


and more useful variation index would be obtained, and one better fitting the 


_ observed distributions. Luigi Gherardellit® j in Italy and Aime Coutagne* in 

France early in the 1930’ s utilized suggestions in their studies of runoff 

_ variations. | Ina brilliant paper classifying streams not only by their indexes 

_of variation with respect to daily stream flow as in the present paper, but also i 

_by a series of other variation indexes to take account of seasonal and subsoil | at 

"seepage characteristics, Mr. Coutagne presents values for the Lane-Lei variation 

coefficient (in terms of its reciprocal) for a number of European streams. 
--values are very similar to those in the paper but his conclusions are different. 


Discussion by John S. Cotton of Reservoirs: Symposium,” ASCE, 
“Statistics by Intercomparison, with Remarks. the Law of Frequency of by 
Galton, The London Edinburgh and o— Piadesophiods Magazine and Journal of Science, 4th Series, Vol. 
XLIX, January to June, 1875, p. 33. 
Sujet de Quelques Statistiques les Déterminations Hydro- 
7 logiques,’”’ by Luigi Gherardelli, Rapports Bulletin No. 21, Assn. Internationale d’ Hydrologie Scientifique, L- 
Rome, 
3 47 “Classification i Cours D’Eau D’Aprés Les Méthodes Statistiques,” by Aime Coutagne, Pithiviers, 7 
Imprimerie des Caisses d'Epargne, Grenoble, France,1935. 
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STREAM 
_ Presumably Mr. Coutagne utilized endent: methods similar to those presented _ 
by the American Society of Testing Materials in computing t the variation 
index, procedures which the | w riter would recommend. Tt is s hoped, i in their 


example, a as it is important that, their method be fully understood those 


their w work and the present procedures are not clear. 


should help ite use. se. A stream “whose logarithmic variability i is 0. 
(the typical stream of Table 2) h has an antilogarithm of AL The theory 
of statistics indicates that a stream with this standard deviation of its” 
~ logarithms « of flow has two thirds of its daily his wilde the range of one 
- fourth to four times the yearly daily average. . Average flow here means the 


= corresponding to the _antilogarithm of the average logarithm, a value of 


the order of 0.6 times the ordinary arithmetic m mean flow. _ Such a stream for 

— 83% of the time would develop 25% of its potential power without storage 

The wide 1 range in the suitability of the streams for ‘power developments from . 

_ purely hydrological viewpoint is indicated by the wide range | of the authors’ 

indexes, from 1.2 to 0. 2. These, having antilogarithms of 16 and 1.6, mean - 
that for 83% of the time these streams had flows that were ‘greater. than n 6% 


and | 63% of the yearly average. 


‘The second part o of the ‘paper | concerns a discussion of t of the factors on whic 


4 cover) and presence and swamps are the = 
_ important factors in stream flow variability during the greater part of the 


time and * * * other factors have relatively smaller effect.” __ 
In view of the n many varying characteristics included in the watersheds of the 
stream cited in the paper (such as rainfall, ‘drainage area, geography, an 
culture—all « of which are commonly thought to effect the uniformity of stream 
flow), it is regretted that the paper does not contain supporting analysis. for 
this conclusion. In Table 2 the authors propose that the variation index lor 
an unfamiliar stream be estimated by utilizing the average variation index of 
— 0.60 typical of their data and then modifying it to values of from 0.7 to 0.2 (a 
ratio of 3.5 to 1) purely on the basis of effects due to varying geological { factors—_ 
that i is, percentages of glacial moraine and bed rock, or varying percentages of | 
lakes or swamps. Other factors are given ni no consideration. writer 
_ proposes briefly to examine whether other factors may not be of equal i impor- 
Rainfall. —That varying rainfall cau causes a is s almost. too well 
_ aie n to. require comment. = Correlations greater than 0.90 are commonly 


found when comparing» these variables whether on the basis of daily, 


2: 


= 


- weekly, monthly, or annual flows. If this is true : a a glance at typical rainfall 
iy variation indexes should give a # guide to the ve variation to sn expected of von 


eser ntation of Data,” Committee E-1, A.S.T. M., Philadelphia, Pa. Pa. April 19 a al 
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for certain European stations: : 


"Lyon, France. 


France. 


=. 


runoffs varying ina manner. ‘The writer suggests that, as more 
daily rainfall indexes become available from locations in the tropics to the — 
thunderstorm belts: of the arid west, the range of rainfall variation will be 


found quite comparable to the range in stream flow va variation presented 


Altitude. varying altitude « affects the runoff 


Contours of 

equal variation 

in annual rainfall 

(Coefficients of 
variation) 


as 
of 


went from the plains to. o the mountains and ascribed it to the more uniform 
& om mountain rainfall. oe ‘That this effect exists: is also sh shown in Fig. 13, from data 
 suj supplied by Winters P. Nakamura.® Here i in terms of the arithmetic deviation | 
to which, roughly, the authors’ index is proportional, the annual rainfall 
variation index changes from 50 to 20, a ratio of 2.5 to 1 as the altitude becomes 


ad Area.—The effect of the area of watershed on the uniformity of its drainage ~ 


: 4 is likewise common knowledge. Big rivers tend to be much less variable than 


49 ‘Rainfall Coefficient of Variation,’’ by Winters P. December 
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| AM FLOW 
Rainfall a ron is 1.¢ as variable as at an Alps station, and, pre- 
y 
_ they found that their index varied from 0.62 to 0.25 (a ratio of 2.5 to 1) as they — | 
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sq miles, the variation index is reduced on the average from 0.76 to 0. ~~ 
to 2 to 1. data, sk shown as follows, indicate similar 
(eq miles) Coefficient 
38000 
2,500 


Rhone, Lyon, France... .. 
Beaucaire, France. . 0.30 
0.29 

‘Nile, Egypt (compiled by 


To conclude, the foregoing data, limited as they are suggest that varia- 4 
ro in rainfall, altitude, or area may each have effects on the “flow variability 
index” comparable to that ascribed by the authors to ‘ “geology. The 
- feels that, before the authors’ conclusion that “ ‘geology”” i is & principal factor in 
runoff variation can be accepted, correlation studies must be made to se 
the effects of variable rainfall, elevation, and area as just pointed out. © 
though some engineers, lacking other data, may wish to accept the aoa : 


recommendations and estimate unknown stream flow variations following the 
racine recommendations ¢ of Table 2, the writer believes that, for a an . unknown 


- to the rainfall variability index, measured or estimated, or if, on the basis of © 
observed physical of a drainage basis, a variation index 


E. McKezg, M. ASCE —It i is gratifying to note that Messrs. Lane and 
Lei have applied s standard statistical procedures to stream flow data and have 
thus analyzed all reliable records of streams east of the 100° meridian. ‘h Too 
often, civil engineers are prone to develop methods and procedures for analyzing 
data with utter disregard for established statistical parameters which are quite 
applicable; and as a result they are plagued with hybrid indexes, coefficients, 
averages, and ratios. “ ‘The authors have attempted wisely, to avoid such hybrid 
parameters by using the well-established probability analysis; | but (as will be 
shown) they have apparently overlooked a significant feature of their “vari- 
ability 
- By its definition and by its — of computation, the variability index 
& be a logarithm. iti is defined by the authors as the standard a 


of the logarithms. of the discharge values (X) at 10% intervals, as taken from the 


- duration curve. As: such, it is nothing more than the logarithm of the geo- 

lop = X — M 


(5° Associate I Prof. of of San. Eng., California Inst. of T Technology, Calif. 


small flash streams. The authors’ data, when plotted against drainage area, 
na ate tha he ao he wate hed in pase NO an miloac to 10 
‘ha 
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i in _—e Xi is a 2 value read from the duration on curve; M, is s the geometric mean, - 
since log M, is the arithmetic mean of the log X values; N is the number of 
observations (in this case 10); and o, is the geometric standard deviation. In 
the form of its logarithm, the geometric standard deviation has little or no 
7 actual s significance, other than as another index to be carried in the minds of 
: engineers. a As the antilogarithm. of the e variability index, however, the g geo- 
- metric standard deviation proper | has a definite meaning and can be applied © 
directly to stream flow data. 
Inasmuch as the duration curves of stream flow data that plot as straight | 
_ lines on on logarithmic probability paper, between limits of 5%} and 95%, may be 
considered. as approximately geometrically ‘normal, the geometric standard 
deviation, o,, is the accepted statistical measure of dispersion (or variability) . 
of such date. _ It has a physical significance in that, when multiplied by, or 
"divided into, the geometric mean (or 50% value in this case), it results in direct | 


3 - values having the same e dimensions as the mean. The product and quotient so 7 in 


obtained represent the values at 84. 13% and 15.87%, respectively. Thus, if 
the geometric mean of certain runoff data is expressed as 1.20 cu ft per sec per F 
sq mile and a, is found to be 2. 20, then the runoff may be expected to be less a 
than 1.20 X 2.20 or 2.64 cu ft p per sec per sq mile for sia 13% of aad eee and less ~ 


than or 0.545 cu ft sec per sq_ mile for 1 15. of f the time. 


4 variability used the authors, which be log To, OF 0.34 in this 


example, has no such ready application. ; Few human minds are accustomed 


to to thinking logarithmically, and, hence, few engineers can envision what 0.34 _ 


means in this instance. 
Although the variability index i is definitely as a by 
its definition and by its magnitude i in the various tables and examples, and ; 
although the characteristic values of the variability coefficient given in Table 2 ’ 
are logarithms, the au authors were apparently thinking of the antilogarithm of 
‘their index, or the true geometric standard when stated (under 


heading, “Use of the Variability Index”): 
“The shape of the duration curve can be obtained by denen: a . straight-— 
— line duration curve on logarithmic probability paper with a slope such that — 
the ratio of the discharge exceeded 15.87% of the time to the discharge 


exceeded 50% of the time was equal to the variability index selected.”’ 
The quoted statement does not hold for the logarithms used in this paper as the 


variability index. It holds only for the true geometric standard deviation. 
: — The writer concurs with the authors in their reasoning that a measure of 
variability of generalized stream flow data is needed and he feels indebted “a 
_ them for the mass of data that has been thus compiled, analyzed, and compared. a 
‘He recommends, however, that the variability index be called by its established — 
statistical name—the e logarithm | of the geometric | standard | deviation—and th that. 3 
this variability be expressed directly rather than as a logarithm. 


Mz. ASCE.— —The 


authors a are e to be on attempt a stream flow 


Design Engr., Safe Harbor Water Power Corp., , Baltimore, Md. 
wtipdeegngtis Engr., Pennsylvania Water & Power Co., Conestoga, Pa. 
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AN ON STREAM FLOW 29 
index. accurate of flow variability, by some some kind 
of rational ‘numerical coefficient, has long been needed i in hydrologic investiga- as 
any step in that direction is salutary. ‘Further commendation is. 


warranted for the authors’ astuteness in utilizing, as a . basis for their mathe-_ 


pal 
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"PERCENTAGE, OF TOTAL TIME FLOW IS EQUAL TO OR LESS THAN 


a 


Cc. Krumbein, F. J. Pettijohn, and G. H. Otto, a as s applied i in 1 the field of sediment _— 
‘particle analysis. It is often more difficult to’ visualize the adaptability of 
ye oe ‘in a foreign field to a given problem t than ‘to devise a method of attack — a 


the problem from fundamental concepts of physics and e1 engineering. Simi- 
larity of approach to two pees is strikingly apparent only when the two 
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been reduced to dimensionless partial ‘differential equations, 
a 


and a relationship of of the factors involved i is at an uncomplicated mathematical 
-. However, as in the case of » any new theories, : all. concepts that are in- — 


cluded must be explained. 


_ The three run- -of-river hydroelectric at Safe. Harbor, Pa., | Holtwood, 
Pa, and Conowingo, Md., form an excellent chain of stream. odes stations 
for ‘the Susquehanna River discharge. _ Although the three plants have Tela-— 
i. tively little storage, all flow data have been reduced to natural conditions, and — 


- > plants themselves can then be considered as having no effect in equalizing _ 


_ stream discharge and fluctuation in flow. Independent flow records at wal 
It was anticipated that, with such data, ‘valuable confirmation ‘might 1 be 7 
- offered for the authors’ studies of the 1 river; +s and it is in fact due to this an anticipa- : 
tion that the following is presented. — For the period from 1917 to 1946, the 
— flow variability index for the Susquehanna River at Holtwood is 0.432. This © - 
value differs but. slightly from the value of 0.44 determined by Messrs. Jae 
and Lei for Harrisburg, Pa., and may easily be accounted for by the fact that - 
 Holtwood i is 47.0 miles downstream from Harrisburg and has a drainage 
3 700 sq miles greater (26, 800 sq miles as compared to 24, 100 sq miles). 7 
Moreover, the geology of the intervening territory is of highly different char- 
acter” from that of most of the drainage basin, “since at Harrisburg the 
Per precominantly mountainous terrain and enters rolling farm country. ot In 


addition to the yearly duration curve, , monthly duration curves for the same 
7 _ period have been prepared and some are depicted in Fig. 14. Using the method - 
by the authors, fi flow variability indexes for ‘each h month were devel- 


oped from these duration curves. results, with aver average ‘monthly flows at 
 Holtw rood, are as follows: 


— 


Flow viability ver: ow 

0306 | 83, 300 

0279 18,600 
.. 0.3824 12 =. 

18,800 


ere indexes vary from 0.265 for ‘April | to 0.415 for oe they present 


an unusual condition, although entirely possible in in this wherein the 
annual variability index seemingly bears no relation to the monthly indexes. i 


Since it can logically be assumed that the geological structure, soil, — a 


: raphy, a and the arrangement of tributaries, with r re 


gard to time of concentra- 
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on tion of surface flow, are cons nstant for any particular river, the eutaiess in the 
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monthly indexes must arise from one or more other variables. 7 An effort was 

made to see if such variation could be accounted for by some function of the - 

weather. The “monthly variability indexes are in Fig. 15 showing 


average temperature, precipitation, and runoff. The | temperature and 
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15 5.—Montay PRECIPITATION, RuNoFF, TEMPERATURE, AND 
Invex For THE Susquewanna River Asove Hottwoop, Pa, 


cipitation as, are based on records” from 33 stations and 58 stations, 
respectively. It i is strikingly evident that no direct correlation between any 7 
of these quantities s and the flow variability exists; nor is my correlation between — oe 
_ the variability index and the runoff coefficient apparent. — _ Temperature appears } 

to be: an important factor and vegetation density would roughly 


The wide variations in the shape of the duration curves, and, hemes, the 
variability indexes, probably stem largely from the effects of snow and frozen . 
ground in delaying runoff during» the winter months: and from the 
percentages of flow contributed by ground | water during the late summer and 
fall months. . Although the deviations in the monthly variability indexes 
neither prove nor disprove the authors’ contention that geology and topography - 
7 are the principal factors involved, they act as a warning against the indis- po * 
eriminate use of indexes to develop duration curves. The foregoing comparison 
a the telatively constant indexes for April, May, and June, with the widely 
"different duration | curves for these months in Fig. 14 illustrates the degree of 
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ot — been mentioned by “Messr s. Lane a and Lei that the topography and © 
a geology | of Pennsylvania are e highly diversified, and that extremely varied flow 
indexes might be expected on individual streams. it should be pointed out 
that, solely because of this variance, any attempt to predict a variability factor b> 
for a given stream would meet with difficulty. As an example of the influence 
; of geology, the authors cite the case of two adjacent small Pennsylvania streams, _* 
Muddy Creek with an index of 0.29 and Conewago C. Creek with a value of 0. 58. 
As the authors mention, Conew: ago Creek drains an area of Triassic sandstone a 
with traprock intrusions, whereas the drainage basin of Muddy Creek is ina 
_ highly laminated and inclined absorbent albite chlorite schist which, were it 
4 not for the laminations, would indicate a high index for Muddy Creek. Such > 
- unapparent geological al considerations as | laminations s would | be completely over- 
- looked in any but the most extensive investigations of a drainage area, and it _ 
- is very likely that, if such comprehensive geologic data were at hand, adequate 
stream 1 records also would be available. 
examination would for, if is an indication 
of the underlying rock structure, identical index values would be predicted for 
the two streams. The Pennsylvania Gazetteer of Streams® lists both streams. 
; as being i in rolling agricultural country with a a See valley flanked by 


of the critical is mentioned by the authors | themselves— 
that is, the case of Mad River in Ohio. ° _ It should be noted also that the period i: 
= “a recorded flow observations for both -Conewago Creek and Muddy Creek on 
only 10 years s duration each; and in addition, in both the 


it is hoped that the authors and others will continue in their efforts to place 
& stream flow variation more nearly on an absolute and mathematical | basis. 
wW. LANE, M. _ ASCE is not possible for Mr. Lei to join the writer 
a preparing the closing discussion of this ] paper because of the disrupted | com- — 
munication betw een China and the United States. _ Although it is Langit oat = 


effect « of ground water, ‘and other factors peculiar to a a given 


_ The large number of constructive criticisms presented by the various dis- 
cussers has added greatly to the value of the paper and they ar are » appreciated. 
_ They fall into two general classes: (1) Those dealing with the statistical aspects 
of the paper and the mathematics of probability paper, and (2) those dealing 
with the hydrological aspects and applications of the variability indexes. ce ww 
- Considering f first the discussions i in class (1), it is advisable to give the reasons 
for adopting the method of handling the flow of data used in the p paper. . The ( 
data used were drawn from a large number of sources—and were given in four 
_ different forms: (1) Plotted duration curves; (2) discharges exceeded during 


88 “Water Resources Inventory Report: Gazetteer of Streams, Pt. III,’’ Water Supply Comm. of 


“ Hydr. Consultant on Chf. Engr.’s Staff, Bureau of Reclamation, Denver, Colo. 
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LANE ON STREAM FLOW) 
7... certain selected percentages of time; (3) percentages ges of time .e that ; the « discharge _ 
exceeds certain selected vile: (4) mean daily ‘discharges: for each day 
covered by the record. | Tous use all these data, a method which would uld handle data 7 = 
expressed i in any of these forms had to be selected. Since future u users of the | - 
suggested variability index will also encounter data i in all these forms, in order a ae 
that they | may all use the same procedure, it is necessary that the method } ; 
selected be ap applicable to all types of data. 7 _ The method of preparing a duration — a 
curve, and of reading vs alues from it at certain selected typical points, | was the 
only | process the writer and Mr. Lei could devise w hich met these requirements. 7 
In making duration-of-flow studies from records | ‘of mean daily flow, the | 
7 ‘shortest m method s seems s to be to divide the entire range of of the stream flow a - 


7 tallying the number of days f found i in each partial range. ” From the numbers sy 


= 


shown by the tallies of these partial r ranges the number of days with sing Be 


below the upper limit of each partial range can be easily obtained by adding the 
_ tallies i in that range to those in all partial Tanges cov ering lower r discharges. : 


‘ee curves, where the data used were in the form of mean daily aaa 
they believe that any other method i is “unnecessarily laborious. 
Mr. Ownbey expressed the opinion that in computing the index the sum of _ 
the squares of the differences of the logarithms fr from the mean ofthelogarithms 
a should be divided by 10 instead of | by 9, as was ; done in the paper. % — 


selecting 9, the writer and Mr. Lei consulted several people who are well 


versed in statistical methods, in addition to studying the methods of handling 
similar data recommended by W . C. Krumbein. _ Still others who are well _ 
acquainted with ‘statistical | theory have discussed this paper without t raising 
‘this question; therefore, it seems likely that the use of 9 by the writer and oe 
Lei is sound. If it should not be sound, however, it w would have no effect on 
the practical usefulness of the index developed in the paper, w hich differs by 
a constant ratio from the one which Mr. Ownbey would obtain. > Say 
‘ow Mr. McKee suggests that the geometric standard deviation would be a 


_ better v value to use for the Vv variability index than the value adopted i in the | paper. 7 


‘the advantage of more re easy designation, on, and, therefore, might be 1 more re accept- 7 
able to statisticians, it would have no practical advantages over the one in the 

paper, which is the antilog of the one Mr. McKee suggests. 
The . use of the coefficient of variation, and the use of discharges e expressed as _ 
ratios to the mean flow have been suggested as possible improvements on the 


index offered in the paper. One of the advantages of the index used by _ ; 
writer and Mr. Lei is that, by using the logarithms, the index obtained gives an 
_ expression of tl the relative variability from the mean, and, therefore, accom- _ 
plishes” the same purpose as the coefficient of variation. Converting the dis- 
_ charges to values expressing their magnitudes in terms of the mean and then de- 
7 termining the variability index would result in exactly the same value as was 
obtained by the method given in the paper. — One of the great advantages of 


ae the index used by the writer and Mr. Lei is that it is independent of the dis- ss. 


units as Mr. Foster has shown. values obtained ar are the 
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‘LANE ON STREAM REAM FLOW 
The w agree with Mr. -Chow’s contention that the variability 
int is a function of the length of the record, when the the time is expressed as a 
"percentage of the total time of flow. = Of « course, there are differences for the 


duration curves of different years, but this is a chance and not a functional ~ 


o an error in the f paper (under the _ 
‘wenden: “Use of the Variability Index”) which resulted from the omission of © 


the word ‘ ‘antilogarithm.’ The correct statement is: 
“The shape « of the duration curve can be obtained by drawing | a straight-_ 
Tine duration curve on logarithmic probability paper with a slope such that 
the ratio of the discharge exceeded 15.87% of the time to the discharge ex- — 
a 50% of the time 1 is equal to the antilogarithm of the variability index 
selected.” 
The mated given by Mr. Foster for obtaining the predicted flow duration é 
£ curve from the estimated value of the variability in index, and the curve e showing 
the ratio of the median to the mean discharge for various values of the varia- 


_ bility i index, as given by Mr. Foster and als also Messrs. s. Ospina and Tama, w ill be : 


of considerable assistance in n applying the methods presented i in the ; paper. The 

_ explanation of the method of constructing logarithmic probability paper offered 

~ by Mr. Foster will be found of great assistance to those who may have coommon 

to make uy up ‘such forms. 
: Turning now to the > hydrological aspects and applications of the variability ; 

*. index, Mr. Wing’s s discussion of the history of the idea of this index is of great 
interest. oa As he has stated, its roots go back over seven decades to Francis 

_ Galton, and the basic idea was s developed by Aime > Coutagne about a decade 
= the independent development of it by the writer and Mr. Lei. | Mr. 


-Coutagne expressed his idea in the form of a uniformity index—the higher hi his 


value the more uniform the flow—whereas the writer and Mr. Lei developed 
8 variability index, the reciprocal of that developed by Mr. Coutagne, with 
* higher values allocated t to the more variable stream flows. rae 
Some e difference of o opinion developed i in a the discussions regarding the varia- 
tion of the index values with the drainage area. - Mr. Lull presents data indi- | 
. very little decrease i in the variability with drainage area; and Mr. Wing | 
and Messrs. Ospina and Tama believe that the variability decreases appre- 
ciably with drainage area. > view ‘that! the? index decreases with ‘drainage 
7 area was | given in the paper. E Howe ever, since ‘the decrease is of a low order 
(a decrease of perhaps about 50% i in index value for an increase of i 
area of 10 0007), it might be argued that the increase was negligible for many — 
It is regretted that the space available did not permit compe 
j just tification of the | contention that g geology and surface storage were the 
most important factors ‘in controlling flow variability within the area cov- 
ered by the data submitted. Had the maps showing the variability indexes _ 


and their relations to local geology and the discussions of watershed conditions — - 
which were a part of the comprehensive paper been ‘published in full, it is Za 


believed that this point would have been much more clearly demonstrated. 
Other factors that been to for the variation were rain- 
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‘LANE ON STREAM FLOW 
fall, variation and Inn regions s of high variations within 
short distances can be appreciable, but it is difficult to believe that such can be _ 
the case in flat States such a as Iowa, , Wisconsin, Illinois, and Ohio In these 


common, and i in n many ¥ watersheds with which this writer ¥ was personally ac- a a 

quainted, neither rainfall nor forests seemed to be p possible explanations. — Mr. . 

Lull gives data which indicate an index about 25% higher for nonforested 


_ watersheds than for forested ones in Iowa and Wisconsin. In these states 


forested areas are allowed to grow timber because they are unsuited to cultiva- 


_ tion, either because of high relief, unsuitable soil, or swamps. . This is the case — 


oh the | Iowa streams” mentioned, and probably also for those in Wisconsin. - 


These conditions favor low variability f for geological reasons, and, therefore, 

‘the lower indexes for r the forested ones care not necessarily due to the forests. 

Perhaps the best arg arguments for the preponderant influence of surface storage 
a the very low values of the St. Lawrence, St. Clair, and St. Marys rivers of 
the Great Lakes “system, which are much lower than any others obtained. 

_ Strong arguments for geology being very important in the area covered by the © 
records are the low flows | of the Kansas and Nebraska streams. _ Although no 
quantitative data on the rainfall variability i in these Great Plains states are 4 

available, there is every reason to believe that it is less uniform than in the 
states farther east. _ The watersheds of these streams are nearly treeless. r a 
spite of both these factors, these Great Plains streams give low indexes. Bo ne 
it It should have been emphasized more forcefully perhaps that the conclu- 
sions of the paper were confined to the areas covered by the flow w records avail- 
o able, and that, until data from other regions are available, it is not possible to 
determine how applicable they are elsewhere. With the exception of some 
streams in North Carolina and perhaps i in southern Illinois, , No cases we were re found | 
which seemed to indicate a rainfall variability effect. However, this writer — 


ag does not claim that rainfall effects are everywhere generally negligible, a 
z only that they were generally negligible in the area covered by the records. In 


1 fact, where there are great differences in rainfall variability there are certain to 
‘be differences in flow variability. The flow variability is the e result of the 

| interaction of rainfall variability and natural storage, either above ground or > 
below ground; the flow variability varies in the same direction as the rainfall 
variability and in an inverse direction from the natural storage. It is Leeper 


Messrs. ‘Bruin and | Hudson have mentioned the necessity of expressing 


that the low values obtained by Messrs. Ospina : and Tama for the streams of — 


total l elapsed time. In all the flows used by the writer and Mr. Lei , these v were | 
the ss same; but in n small streams, where the flow ceases a part of the time, it is. 


— 


necessary to use only the time of flow. 
pi Mr. Kazmann has called : attention to the lack of exactness i ‘in the terms arin 


‘in 
_ exactness would be desirable, but the labor involved in working up the data on 


individual watersheds, in quantitative t terms on pom ch generalizations could 


< 
= 


‘= 
| 
— 
— 
| 
j j 
1 
| 


= 


— vance would h have been prohibitive, even if the data | had been available. — In 
cases this was not true. 
serve as a in handling this variable. 


SOME LAKE. AREAS 


Messrs. Schroyer -and Schuleen have shown the great differences in varia-— 


- 4 bility indexes of flow for the Susquehanna River in the different months of the 


4 year. 7 This serves to emphasize the necessity of having the data cover one or 
‘more ‘complete cycles, such as years, in which the conditions are much the same, = 


rather than to compare conditions which are consistently different, such as 
those i in one month of the year as compared to ) those in another month. ; 


In judging the geological and surface characteristics influencing flow 
variability in drainage areas for which geological data are not available, 
photographs. or airplane reconnaissance would be very valuable. Considerable a 
work been toward interpreting these conditions from air photo- 


‘graphs, 66,56,57,58 which would repay study by those concerned. 
‘‘Aerial Photographs Their by Harold T. U. Smith, D. Appleton- Century Co., 
Ine., New York, N. Y.,1948. 
$7 Origa, Distribution, and Airphoto Identification of United States Soils,” by D. S. Jenkins, 
BRS Belcher, L. E. Gregg, and K. B. Woods, Technical Development Report No. 62, Civil Aeronautics Ad- 
ministration, Washington, D.C.,1946. 
- m & ‘‘Aerial Photographs Used for an Engineering Evaluation of Soil Materials,” by R. E. Frost and K. 


Proceedings, 2d International on a Boil Mechanics and Foundation Eng., Rotterdam, 
Vol. I , 1948, pp. 324-330, 
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